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Roof-spray cooling systems have been developed and implemented to 
reduce the heat gain through roofs so that conventional cooling systems can be 
reduced in size or eliminated. Currently, roof-spray systems are achieving 
greater effectiveness due to the availability of direct digital controls. The 
objective of this thesis was to develop a mathematical model of the heat 
transfer through a roof-spray cooled roof that predicted heat transfer based 
on existing weather data and roof heat transfer characteristics as described by 
the Transfer Function Method. The predicted results of this model were 
compared to the results of existing experimental data from previously 
conducted roof-spray cooling experiments. 
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The mathematical model is based on energy balances at the exterior and 
interior surfaces of the roof construction that include evaporative, convective, 
radiative, and conductive heat transfer mechanisms. The Transfer Function 
Method is used to relate the energy balances at the two surfaces that differ in 
amplitude and phase due to the thermal resistance and thermal capacitance 
characteristics of the roof. 

The combined model yielded moderately good predictions of heat 
transfer through the roof experimental results for the roof-spray cooled 
condition. The calculation method shows promise as a relatively simple means 
of predicting heat gains based on calculation procedures that are similar to 
those frequently used by many engineers. 
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CHAPTER 1 
INTRODUCTION 



Roof-spray cooling systems continue to increase in popularity for a 
number of reasons. Evaporative roof-spray cooling which has been widely 
implemented since the 1930s to provide cost effective cooling for industrial 
applications is now seen as a method to provide cooling which does not add to 
the depletion of fossil fuels or contribute to possible global warming problems 
caused by combustion and the use of CFC based refrigerants in traditional 
vapor compression refrigeration systems. A major advance that has improved 
the effectiveness of roof-spray cooling is the availability of direct digital 
controls, which allow water to be applied to roofs to achieve maximum 
cooling without the application of excessive water. 

Many researchers have studied various forms of rooftop evaporative 
cooling including solar ponds, trickling water, and roof-sprays. In 1940 
Houghton et al. [1] studied the effects of roof ponds and roof-sprays on 
temperature and heat transfer in various roof constructions. In the 1960s, 
Yellott [2] reported on the effectiveness of a roof-spray cooling system that 
consisted of little more than a rooftop grid of sprinkler heads and supply pipes 
with a water supply controlled by a solenoid valve and timer. More recently. 
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research has been conducted to numerically model roof-spray cooling systems. 
This includes the work of Tiwari et al. [3], Carrasco et al. [4], Somasundaram 
and Carrasco [5], and Kondepudi [6], 

The heat transfer through roofs with various construction 
configurations has also been studied extensively. Spitler and McQuiston [7] 
discussed the latest modeling techniques for the prediction of cooling loads 
and times as they occur in the zones of a building as a result of the daily loads 
of ambient temperature and incident solar radiation. 

However, a model that incorporates the current modeling of roof 
systems and roof-spray cooling is lacking. This thesis presents a numerical 
model for the combined heat transfer response to roof-spray cooling and roof 
construction in response to the variation of ambient temperature, humidity and 
solar flux over a typical day and incorporates the effects of thermal 
capacitance of a roof as modeled by the Transfer Function Method (TFM). The 
TFM uses conduction transfer functions to predict the heat flux at the inside 
of the roof as a function of previous values of outside and inside temperatures 
and heat fluxes. The magnitude and direction of the inside heat flux may be 
smaller or larger than the heat flux at the outside surface depending on 
whether the thermal mass of the roof is increasing or decreasing in stored 
energy. This thesis will present a numerical solution to the heat transfer 
problem and compare the results to existing experimental data from previously 
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conducted roof-spray cooling experiments. The heat transfer mechanisms 
acting on the outside surface in response to ambient temperature and incident 
solar radiation are evaporation, convection, radiation, and conduction. The 
transfer mechanisms acting on the inside surface are convection, radiation, and 
conduction in response to the energy transfer at the outside surface and a 
fixed room ambient temperature. The energy balances at the two surfaces are 
related, but the amplitude and phase of the inside energy transfer will depend 
on the thermal resistance and thermal capacitance characteristics of the roof. 

The solution presented should enhance the understanding of the 
combined effects of evaporative cooling and thermal capacitance of building 
envelope construction. As a design tool the benefit of this modeling will be 
that building designers would more easily choose the best combination of 
roofing, spray-cooling, and conventional cooling systems. This modeling also 
has an added benefit in that it can be used as a basis for algorithms that 
determine the most effective and efficient application of roof-spray cooling to 
a roof. 
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CHAPTER 2 

MATHEMATICAL MODEL 



Introduction 



The various elements of the heat transfer mechanisms are presented in 
this chapter. First, the mathematical basis (derived in terms of the energy 
balances) is presented in the form of equations. The different heat transfer 
mechanisms and assumptions concerning these mechanisms are also discussed. 
Second, the thermal capacitance characteristics and effects of the roof system 
are presented. The different types of roof construction, the heat transfer 
properties of the materials, and the assumptions concerning these properties 
are also discussed. Finally the interaction of the heat transfer mechanisms and 
the construction are discussed and a mathematical basis in the form of 
equations is presented. 

Heat Transfer Mechanisms 



Kondepudi [6] and Carrasco et al. [4] presented a thorough 
development of a mathematical model of the heat transfer mechanisms of a 
spray cooled roof with some assumptions being made. These assumptions 
include: 
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(1) The problem is quasi-steady in response to varying ambient 
temperatures and solar flux. 

(2) An infinitesimally thick film of water is maintained on the roof 
surface. 

(3) The saturated water vapor pressures respond in a linear fashion 
to temperature. 

(4) Wind effects are included in the convective heat transfer model. 

(5) Thermal capacitance effects of the roof and water film are 
neglected. 

All of these assumptions with the important exception of the thermal 
capacitance effects of the roof will be carried through this paper. The initial 
formulation of the energy balances will ignore the thermal capacitance effects 
of the roof. 

The energy balances for the two interfaces of the process are presented 
in Figure 2.1. The energy balances at each interface are described below and 
provide a problem that must be solved simultaneously. 

Roof Exterior-Air Interface 

As indicated in Figure 2.1 the incident solar energy is dissipated by a 
combination of evaporative, radiative, convective, and conductive heat 
transfer. 



^soUr ^cv»p"^^r4 



r»d(outsidc)"^^conv(outside) ' 'icond(roof) 






(2.1) 
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Figure 2.1; Energy balances at roof surfaces 



The evaporative heat transfer is a result of Pick's Law of diffusion of the water 
vapor at the wetted surface diffusing into the ambient air, as presented by 
Tiwari et al. [3], 



q =0.013h (P -P ) 

Mev*p ev*pV^ wr w*/ 



( 2 . 2 ) 



The partial pressure of water at the roof surface, P^^, is predicted from 
Mathur's correlation [8] for the water saturation pressure, P^^, as a function of 
temperature. 



+ A,T + + A3T' + AJ“ 



where Aq = 670.2012 psia 
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=-5.325521 psia/'R 
A 2 = 1.59464 XlO-2 psiarR' 

A 3 = -0.2134061 XIO--* psiarR' 

A, = 0.1077853 XIO'"^ psia/^R" 

= saturation pressure of water vapor, psia 
T = absolute temperature, ®R 

The partial pressure of water in the ambient air is predicted by ASHRAE [9] 
as: 

Pw. = (Pco)/(0.62198 + oj) 



where P = total barometric pressure of moist air. Pa 
(j) = humidity ratio of the moist air, 

9rid(oui»ide)“^^(^r,o“Ta) (2.3) 

The equation for the outside radiative heat flux can be refined based on 
Alford et al. [10], who proposed that the ambient temperature for the radiative 
heat transfer should be decreased by 2.8 K. This decrease is based on the 
assumption that all radiation reaching the earth originates in the first 1500 
meters of elevation and accounts for a normal atmospheric temperature 
gradient of about 7 K per 1000 meters elevation. Thus, 



q,.d„u«u,re<^(T?.„-(Ta-2.8)^) (2.4) 

Another method of correcting for the difference in net longwave 
radiation exchange between the roof and the sky is presented by Bliss [11]. 
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Under clear sky conditions the sky-to-surface radiation difference for a 
horizontal black surface, 5R, is typically in the range of 63 to 95 W/m^ This 

radiation difference is a function of elevation, cloudiness, air temperature near 

the ground, and the dew point, but will be assumed to be an average of 63 

W/m^ for clear sky conditions at the experiment sites of Bryan, Texas and 

Pittsburgh, Pennsylvania. Bliss' model will be used in the predicted heat flux 

calculations in this report. 

The convective and conductive heat fluxes at the outside roof surface 
are given by the following two equations: 



Roof Inside Surface-Room Air Interface 

As indicated in Figure 2.1, the heat transfer mechanisms acting at the 
inside surface are conduction, radiation, and convection. 




(2.5) 



9cond(rooO 



(2.6) 




(2.7) 



9rad(insidc) 



( 2 . 8 ) 




(2.9) 
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On the basis of weather data from a database that includes solar 
radiation, ambient temperature, and humidity, the above equations can be 
solved iteratively to obtain a predicted roof surface temperature and a 
predicted conductive heat flux into the roof. This is based on Kondepudi's 
assumption [6] that there are no thermal capacitance effects of the roof or 
water film 

Thermal Capacitance Effects of Roof 

The above discussion provides a simple methodology for predicting the 
outside roof temperatures and heat fluxes if the thermal capacitance effects of 
the roof are neglected. In actual applications the thermal mass effects of the 
building envelope can play an important role in decreasing and/or shifting the 
phase of the cooling load on a building. Several load calculation 
methodologies have been developed to account for this thermal mass effect. 
Spitler and McQuiston [12] presented a review of the most popular cooling 
load calculation methodologies that are being used for manual and computer 
calculations. Because of the proliferation of computers in calculating cooling 
loads. The American Society of Heating, Refrigerating, and Air-Conditioning 
Engineers (ASHRAE) [13] is currently concentrating its research on improving 
and expanding the range of circumstances to which the TFM can be applied. 

The TFM, as described by Mitalas [14], was developed in the late 
1960's as a computer-oriented load calculation method that applies weighting 
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factors to the loads experienced by the exterior surfaces of a building to 
account for the thermal storage effects of the building envelope. These 
weighting factors are referred to as conduction transfer function (CTF) 
coefficients and relate the difference in the sol-air temperature and the inside 
space temperature to predict the difference in magnitude and phase shift 
between the load on the outside surface of the structure and the load that is 
converted to a cooling load. Sol-air temperature [13] is an equivalent outside 
air temperature that, in the absence of all radiation changes, accounts for the 
heat transferred to the exterior surface due to the combination of incident 
solar radiation and radiative and convective heat transfer between the exterior 
surface and the ambient air based on 



T.=T +Qfq.„,„ /h,-e5R/ho 

which is derived from the following two equations [13]; 



( 2 . 10 ) 



+h„(T.-T,.)-£6R 



( 2 . 11 ) 



q=h„(T.-T,„) 



( 2 . 12 ) 



Based on the sol-air temperature described above and a constant inside room 
air temperature, the TFM predicts the heat gain through a roof by the 
following equation. 




f m.. 



1 1 

(2.13) 



CJe.0 S l^n(TeJ-n6) ^ dn(C|e,0-n6) T'set S Cn 

n=D n=l n=0 



The CTF coefficients take into account the heat transfer coefficients at the 
exterior and interior surfaces of the roof and the roof construction. CTF 
coefficients can be calculated using a computer program presented by 
McQuiston and Spitler [7] or use approximate values of CTF coefficients that 
were developed based on investigations by Harris and McQuiston [15] into the 
thermal characteristics of walls and roofs and their effect on the heat transfer 
to conditioned spaces. These approximate values are available in tabular form 
in various texts (see [7] and [13] , for instance) and are included in many 
computerized load calculation programs such as the software provided in 

Carrasco et al. [4] and the Load Design Program [16]. 

Combined Mathematical Model 



The models for evaporative roof-spray cooling by Kondepudi [6] and 
the Transfer Function Method (TFM) described above were combined to 
predict the effect of evaporative roof-spray cooling and the thermal 
capacitance effect of roof construction. Due to the iterative nature of the 
solution of the roof-spray cooling problem and the application of the TFM, the 
use of a computer to solve the combined problem is virtually imperative. The 
problem is solved in a two-step process that first ignores the thermal 
capacitance effect of the roof to predict a heat flux into the roof and then 
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applies the TFM to predict the cooling load on the interior space that results 
from the calculated exterior heat flux. This solution predicts the amplitude 
decay in the load as a result of the evaporative cooling and the roof 
construction and the phase shift in the cooling load as a result of the thermal 
behavior of the roof in response to imposed loads. 

Experimental Roof Parameters 

The thermal characteristics of the roof that Somasundaram et al. [5] 
used in their experiments are provided in Table 2.1, which includes conduction 
transfer coefficients obtained from Reference [16]. The thermal properties of 
each layer were taken from measurements made by Somasundaram et al. [5] 
and from tabular data from Reference [13]. A study of the resulting amplitude 
decay and phase shift in the cooling load is conducted in Chapter 3. 

Weather Model and Data 



The weather data that was used for the prediction of heat transfer in the 
Bryan, Texas roof was taken from the Typical Meteorological Year (TMY) 
readings for Waco, Texas averaged over the period of the experimental 
roof-spray cooling. The data of Waco, Texas were used because it is located 
only 75 miles from Bryan, Texas, which has similar summer design 
temperatures, and it is only 70 meters higher in elevation [13]. The ASHRAE 
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solar heat flux calculation method [13] which is based on the work by 
Threlkeld and Jordan [17] could also be used TMY weather tape was 
selected for its solar radiation data that makes it a good choice for solar 
design problems. 

Actual experiment site weather measurements including ambient wet- 
and dry-bulb temperatures, solar radiation, and wind speed were reported by 
Houghten et al. [1] for the experiments in Pittsburgh. Ambient pressure data 
were extracted from Trane Typical Cooling Day weather files for Pittsburgh 
[16]. Typical Cooling Day pressure varied only 340 Pa, so a constant pressure 
could have been assumed with little effect on predicted results. 
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Table 2.1; Thermal Characteristics of Roof at Bryan, Texas 



Physical Description 



Description 


Thickness Conductivity 
{mm} (W/m^ F} 


Density 

{kg/m^} 


Specific Heat 
{kJ/(kgK)} 


Resistance 
{(m^ K)AV} 


Outside air film 


0 


.000 


0 


.00 


.044 


Rigid insulation 


25 


.043 


32 


.84 


.587 


Felt and asphalt 


10 


.192 


1121 


1.68 


.050 


Gypsum roof deck 


150 


.156 


800 


1.09 


.978 


Mineral insulation 


150 


.069 


10 


.84 


2.201 


Insulation air film 


0 


.000 


0 


.00 


.120 


Decrement factor =120, 


Overall weight =135 kg/m 


1 ^ , Overall U-value = .250 W/(m^ K) 


Time delay = 1 1 Hours, Heat capacity 


= 337 kJ/(m^ kg K), C. 


coefficient = . 


183E-02 



n 


coefficients 


D„ coefficients 


0 


.1733146E-08 


.lOOOOOOE+01 


1 


.1402007E-04 


-.2086841E+01 


2 


.3365204E-03 


.1422010E+01 


3 


.9128446E-03 


-.3563488E+00 


4 


.5027138E-03 


.2923254E-01 


5 


.6378616E-04 


-.7691044E-03 


6 


1774670E-05 


.5973486E-05 
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Table 2.2: Thermal Characteristics of 50 mm Thick Concrete Roof at 

Pittsburgh 



Physical Description 



Description 


Thickness Conductivity 
{mm} {W/m^K} 


Density 

{kg/m^} 


Specific Heat 
(kJ/(kgK)} 


Resistance 
{(m^ K)AV) 


Outside air film 


0 


.000 


0 


.00 


.044 


Slag 


13 


1.436 


880 


1.68 


.009 


Felt and asphalt 


10 


.192 


1120 


1.68 


.050 


Slag 


13 


1.436 


880 


1.68 


.009 


Felt and Asphalt 


19 


.192 


1120 


1.68 


.050 


Concrete 


50 


1.731 


2240 


.84 


.029 


Inside air film 


0 


.000 


0 


.00 


.120 



Decrement factor = .775, Overall weight = 168 kg/m^ , Overall U-value = .346 W/(m^ K) 
Time delay = 4 Hours, Heat capacity = 4 1 0 kJ/(m^ kg K), C„ coefficient - .653 



n coefficients 



coefficients 



0 3732012E-01 

1 .3805050E+01 

2 .2269784E+00 

3 .9101764E-02 

4 .7930927E-05 

5 .8646022E-11 

6 0 



.lOOOOOOE+01 

-.9006944E+00 

.1360358E+00 

-.2234380E-03 

.5431017E-08 

0 

0 
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Table 2.3; Thermal Characteristics of 50 mm Thick Pine Plank Roof at 
Pittsburgh 



Physical Description 



Description 


Thickness Conductivity 


Density 


Specific Heat 


Resistance 




{mm} 


{W/m^ F} 


{kg/m^} 


(kJ/fkgK)} 


{(m^ K)AV} 


Outside air film 


0 


.000 


0 


.00 


.044 


Slag 


13 


1.436 


880 


1.68 


.009 


Felt and asphalt 


10 


.192 


1120 


1.68 


.050 


Slag 


13 


1.436 


880 


1.68 


.009 


Felt and Asphalt 


19 


.192 


1120 


1.68 


.050 


Pine Planking 


50 


.121 


590 


2.51 


419 


Inside air film 


0 


.000 


0 


.00 


.120 


Decrement factor = 


.708, Overall weight = 84 kg/m^ 


, Overall U-value =1.33 W/(m^ K) 


Time delay = 5 Hours, Heat capacity 


= 370 U/(m' kg K), C. 


coefficient = 


189 



n 


coefficients 


D^ coefficients 


0 


.3732012E-01 


.lOOOOOOE+01 


1 


3805050E+01 


-.9006944E+00 


2 


.2269784E+00 


.1360358E+00 


3 


.9101764E-02 


-.2234380E-03 


4 


.7930927E-05 


.5431017E-08 


5 


.8646022E-11 


0 


6 


0 


0 



CHAPTER 3 
RESULTS 



Introduction 



This section presents the numerical results of predicted cooling loads 
along with the experimental results. First, the predicted results for a structure 
in Bryan, Texas are compared to the experimental results obtained by 
Somasundaram and Carrasco [5]. Then predicted cooling loads for two roofs 
at Pittsburgh are compared to experimental results of Houghten et al. [1], 

Roof Heat Transfer and Temperatures at Brvan. Texas 

The comparisons of the experimental data and numerical predictions for 

the roof-spray cooled roof at Bryan, Texas are plotted in Figures 3.1 and 3.2 

and provided in Table 3.1. The predicted average heat gain through the roof, 
q^, of -4.7 ±1.2 W/m^ is in very good agreement with the experimental value 
of -4.2 W/m^. The predicted outside roof surface temperatures are in good 

agreement with the experimental data, with an average overprediction of less 

than one °C. The peak predicted conductive heat flux through the outside of 

the roof, q^onturoof)’ lower than the peak experimental data. The average 

predicted q^ondcrooc 90% higher than the experimental results. Some 
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variations are expected due to the thermal capacitance effects of the roof that 
were neglected in the prediction of the outer conductive heat flux and the fact 
that the outside flux meter was located below the gypsum roof deck in the 
experiments. 

The results plotted in Figures 3.1 and 3.2 also indicate two likely 
problem areas with the models used to predict results. These are weather and 
the thermal properties of the roof. First, although the cumulative daily solar 
radiation obtained from the TMY weather data [17] is very near the 
experimental amount. Table 3.2, the peak experimental solar radiation was 
likely higher at midday and tapered off more rapidly at times other than 
midday. This would account, in part, for the differences in the shapes of the 
predicted and experimental outside surface heat flux and temperature curves. 
The second likely weather problem was an underprediction of the net 
nocturnal radiation exchange between the roof and the atmosphere that was 
reflected in the underprediction of the nighttime cooling effect. This also 
would have contributed to the predicted average q^ondOocf) being 90% higher 
than the experimental average value. Minor inaccuracies in the thermal 
properties of the roof model including the convective heat transfer coefficients 
also would have contributed to some of the differences in the experimental and 
predicted results. This can be illustrated by the fact that the predicted 
difference between the maximum and minimum values of is 0.7 W/m^ while 
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the experimental difference is 1.6 W/m^ which may be a result of inaccuracies 
in the mass and thermal conductivity of the roof materials used in the 
mathematical model of the roof provided in Table 2.1 

As a check on the validity of the roof and exterior load models, 
experimental results from dry, no roof-spray cooling, tests on the same roof 
during the period of 17-21 August were compared to predicted results. The 
comparison of experimental data and numerical predictions for the dry roof are 
plotted in Figures 3.3 and 3.4. Due to the 21% difference in experimental and 
TMY cumulative solar radiation data, the TMY hourly solar radiation data 
were adjusted upward to match the experimental daily cumulative data. The 
TMY hourly ambient temperatures averaged only 1.1 X lower than the 
experimental data and were not adjusted. 

The predicted difference between the maximum and minimum values of 

is 1.2 W/m^ while the experimental difference is 1.3 W/m^ which indicates a 

good representation of the thermal capacitance effects of the roof. The 
predicted average of 2.4 ±1 W/m^ is in moderately good agreement with 
the experimental value of 0.8 W/m^. This difference may be a result of a 

combination of inaccuracies in the weather, heat transfer model, and thermal 

characteristics of the roof that are magnified by the higher heat fluxes 

associated with the dry roof condition. 



ROOF SURFACE TEMPERATURE (DB3 C) 
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Figure 3.1; Relation Between Time and Heat Flux Trough Outside and Inside 
Surfaces of Roof-Spray Cooled Roof at Bryan, Texas, in July 




Figure 3.2: Relation Between Time and Temperature of Outside Surface of 

Roof-Spray Cooled Roof at Bryan, Texas, in July 



SOLAR FLUX ON ROOF SURFACE (W/m^) 






3 
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Table 3.1: Experimental and Predicted Heat Flux and Temperature 

Comparison for Roof-Spray Cooled Roof at Bryan, Texas 



July 17 - July 30, 1987 
Outside surface Inside surface 



Experimental 


Predicted 


Experimental 


Predicted 


Average Roof Heat Flux 
(W/m') 


0.40 


0.76 


0.44 


0.40 

±0.2 


Maximum and Minimum 


7.1 


3.9 


1.3 


0.8 


Roof Heat Flux 
(W/m') 


-3.1 


-2.1 


-0.2 


0.1 


Cumulative Daily Roof 
Heat Flux 
(W h/m') 


9.5 


18.3 


10.5 


10.0 


Average Roof Surface 

Temperature 

CO 


28.2 


29.1 


26.7 


Note; 1 


Maximum and Minimum 


41.4 


41.2 


27.7 


Note; 1 


Roof Surface 


20.3 


18.8 


25.8 


Note; 1 



Temperature 

(“C) 



Note: 1. Experimental data for inside surface temperature were used in the 

predict other table values. 



Source; Experimental data from reference [5]. 
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Table 3.2 Experimental and Predicted Weather for Bryan, Texas, and 
Pittsburgh 





Bryan, Texas 


Pittsburgh 




Wet Roof 


Dry Roof 


Wet Roof 




Test Date 


Test Date 


Test Date 




July 17- 


August 17- 


8 September- 




July 30 


August 21 


9 September 


Experimental Predicted 


Experimental Predicted 


Experimental 


Cumulative Daily 


6398 6497 


6492 5114 


3613 


Solar Radiation 




Note: 1 




(W hr/m^ /24 hrs) 








Average Ambient 


28.54 28.64 


30.24 29.12 


26.3 


Temperature 








(°C) 








Maximum and Minimum 34.79 34.44 


37.49 33.89 


35.0 


Ambient Temperatures 


23.12 22.72 


24.33 25.00 


18.6 


(°C) 








Average Wind Speed 


Note: 2 14.2 


Note: 2 13.6 


22.3 


(k/hr) 








Maximum and Minimum Note: 2 17.8 


Note: 2 16.7 


38.6 


Hourly Average 


Note: 2 10.2 


Note: 2 10.6 


9.7 



Wind Speed 
(k/hr) 



Notes: 1. Due to large difference between predicted and experimental values, 
solar data were adjusted to experimental value for predicting other 
table values. 

2. No experimental wind data available. 

Sources: Experimental data from Reference [5]. 

Predicted data from Reference [17]. 



23 




800 



600 



O 

O 

2 

o 



o 



Figure 3.3: Relation Between Time and Heat Flux Trough Outside and Inside 

Surfaces of Dry Roof at Bryan, Texas, in July 




Figure 3.4; Relation Between Time and Temperature of Outside Surface of 
Dry Roof at Bryan, Texas, in July 
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Table 3.3: Experimental and Predicted Heat Flux and Temperature 

Comparison for Dry Roof at Bryan, Texas 



August 17 - August 21, 1987 
Outside surface Inside surface 

Experimental Predicted Experimental Predicted 



Average Roof Heat Flux 
(W/m') 


3.4 


2.4 


0.8 


2.4 

±1.0 


Maximum and Minimum 


17.6 


9.7 


1.5 


3.0 


Roof Heat Flux 
(W/m^) 


-4.1 


-1.4 


0.2 


1.8 


Cumulative Daily Roof 
Heat Flux 
(W h/m') 


82.5 


58.5 


19.3 


58.2 


Average Roof Surface 
Temperature 

ro 


38.8 


37.2 


28.3 


Note: 1 


Maximum and Minimum 


66.7 


63.6 


29.4 


Note: 1 


Roof Surface 
Temperature (“C) 


21.1 


23.1 


26.9 


Note: 1 



Note: 1. Experimental data for inside surface temperature was used in the 
prediction of other table values. 

Source: Experimental data from reference [5], 



Roof Heat Transfer at Pittsburgh 



The numerical predictions and experimental data for the concrete roof, 
which are given in Table 3.4 and in Figure 3.5, show good agreement. No 
experimental outside roof heat flux or temperature data are available. 
However, the experimental and predicted inside heat flux for the concrete roof 
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are in good agreement in hourly magnitude, phase shift, and cumulative daily 
heat flux into the room. The predicted average of -4.7 ±1.2 W/m^ is in very 
good agreement with the experimental value of -4.2 W/m^ The predicted 

difference between the maximum and minimum values of is 28.2 W/m^ 34% 

higher than the experimental difference, 21.1 W/m^. The differences in the 

predicted and experimental inside heat fluxes may, in part, be attributed to 

minor inaccuracies in the mathematical descriptions of the roof and the models 

of heat transfer at the roof surfaces. 

The numerical predictions and experimental data for the pine plank 

roof, which are given in Table 3.5 and in Figure 3.6, show moderately good 
agreement. The average predicted q^ of -4.2 ±1.4 W/m^ is in moderately good 
agreement with the experimental value of -2.6 W/m^ The difference between 

the predicted and experimental values and the approximate two hour difference 

in the phases of the predicted and experimental q^ indicate that the 

mathematical description of the pine plank roof is inaccurate. The 

mathematical descriptions of the Pittsburgh roofs are based on sketches by 

Houghten et al. [1] that provide no thermal characteristics of the roof 

materials. There is a large uncertainty in the thermal response modeling of the 

pine plank roof because of the unknown characteristics of the pine planking. 

It is unknown if the two inch thickness of the pine planking given by Houghten 

is the actual or nominal thickness. Also, the density, thermal conductivity, 

and specific heat of the wood could vary by over 30% depending on the 
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specific species and moisture content of the planking. An example of the 
variation in the thermal properties of a material from expected is the thermal 
conductivity of the Bryan, Texas mineral insulation that testing [5] revealed 
was 60% higher than values found in Reference [9], 

Summary 

The computer code, based on the mathematical models, provides a 
moderate to good comparison with experimental results of heat flux through 
the inside surface of a roof and into a room. The predicted average heat gain 
through the roof for the Bryan, Texas insulated concrete roof, the Pittsburgh 
concrete roof, and Pittsburgh pine plank roof are in good agreement with the 
experimental averages. The prediction of heat flux and temperature at the 
outside surface of the roof provided good qualitative comparison with 
experimental results. However, large magnitude and small phase differences 
from experimental data were obtained as a result of the thermal capacitance 
effects of the roof that were neglected in the prediction of the outside surface 
temperature and heat flux and the placement of the outside roof heat flux 
meter below the gypsum decking. 
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Figure 3.5; Relation Between Time and Heat Flux Through Outside and 
Inside Surfaces of Concrete Roof at Pittsburgh 
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Figure 3.6; Relation Between Time and Heat Flux Through Outside and 
Inside Surfaces of Pine Plank Roof at Pittsburgh 



28 



Table 3.4: Experimental and Predicted Heat Flux Comparison for 

Roof-Spray Cooled Concrete Roof at Pittsburgh 







August 8, 1939 






Outside surface 


Inside surface 




Predicted 


Experimental 


Predicted 


Average Roof Heat Flux 


1.0 


-4.2 


-4.7 


(W/m') 






±1.2 


Maximum and Minimum 


15.6 


6.3 


9.1 


Roof Heat Flux 
(W/m') 


-12.9 


-14.8 


-19.1 


Cumulative Daily Roof 
Heat Flux 
(W h/m2) 


24.6 


-99.5 


-113.5 


Source: Experimental data from reference [1]. 



Table 3.5: Experimental and Predicted Heat Flux Comparison for 

Roof-Spray Cooled Pine Plank Roof at Pittsburgh 




August 8, 1939 






Outside surface 


Inside surface 




Predicted 


Experimental 


Predicted 


Average Roof Heat Flux 
(W/m') 


0.9 


-2.6 


-4.2 

±1.4 


Maximum and Minimum 


16.4 


2.4 


7.1 


Roof Heat Flux 
(W/m') 


-11.9 


-8.5 


-15.9 


Cumulative Daily Roof 

Heat Flux 

(Wh/m') 


21.2 


-62.3 


-101.7 



Source: Experimental data from reference [1]. 



CHAPTER 4 

CONCLUSIONS AND RECOMMENDATIONS 
Conclusions 

The objective of this thesis was to combine existing models for 
evaporative roof-spray cooling and the transient heat flow through roofs to 
provide a model that will predict roof cooling loads based on known weather 
data and roof material descriptions. The best evaluation of the effectiveness 
of the model is by comparison to well-documented experimental data. 

Two roof conditions were studied. One was dry roof (no spray 
cooling), and the other was damp roof (spray cooled). Three different roof 
constructions at two locations were studied. The roofs included an insulated 
gypsum deck, an uninsulated concrete deck, and an uninsulated pine plank 
deck. 

The mathematical formulation gave a simple but effective prediction of 
the roof cooling load for an interior that is maintained at a fairly constant 
temperature. The agreement of experimental data with predictions of roof 
cooling loads as a result of the combined effects of roof-spray cooling and 
roof construction gave a high level of confidence that the use of the combined 
roof-spray cooling and Transfer Function Method model will yield good 
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predictions for buildings maintained at a constant internal temperature. The 
agreement between predicted and experimental average heat gain through the 
roof, q^, values was very good making the predicted values a useful tool in 
predicting the energy usage for roof-spray cooled buildings. 

The uncertainty analysis described in Appendix B indicate that the 
thermal characteristics of the composite roof are important factors in dry and 
damp roof cooling load calculations and must be accurately described to 
ensure the accuracy of the predicted heat gains through the roof . Two other 
variables that have a large affect on the uncertainty of the predicted heat gain 
are the external convective heat transfer coefficient, h^, and the solar radiation 
transferred to the roof surface as a result of the incident solar radiation and 
the absorptance of the surface of the roof Since h^ is largely a function of the 
wind speed across the roof, accurate'measurement of the wind speed at the 
building location is important to allow the accurate prediction of h„. The 
effect of changes in the absorptance of the roof surface indicates the 
importance of this measurement to allow the accurate prediction of solar 
radiation transferred to the roof. It also indicates the value of low 
absorptance roof surfaces for dry and roof-spray cooled buildings. The need 
to accurately measure the wind speed and incident solar radiation is important, 
not only for experimental activities, but for use in the prediction of the 
required application of water to the roof surface. Therefore the installation of 
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anemometric and heliopyrometric devices as part of the roof-spray cooling 
system may be desirable for the optimization of roof sprays for larger 
buildings. 



Recommendations 



The following is a list of suggested improvements that can be included 
in the modeling of roof-spray cooled structures. 

1) The variation in wind speed over the roof can be incorporated into 
the convective heat transfer portion of the outside load calculation. 

2) The effect of variations in cloud cover and sky clearness on the net 
radiative heat exchange between the roof and the atmosphere can be included 
in the outside load calculation. 

3) The predicted evaporative heat transfer rates can be extended to 
predict the variation in the required rate of delivery of water to the roof to 
optimize the cooling effect without the application of excess water. 

4) Various roof surface materials such as aggregate ballast of different 
sizes and colors can be compared to smooth roof surfaces to determine any 
effects on the heat transfer at the roof surface. 

The mathematical formulation provided acceptable results for buildings 
maintained at constant inside temperatures. The validation of this model for 
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buildings that are maintained at different temperatures during a daily or 
weekly cycle is needed because of the common usage of temperature setbacks 
as an energy conservation measure. 

Any further experiments on roof-spray cooling should include a study of 
the roofs thermal properties including surface absorptivities and emissivities 
in addition to ambient weather conditions that include dry-bulb temperature, 
humidity, wind speed, pressure, solar radiation, cloud cover, and precipitation. 



APPENDIX A 

COMPUTER CODE LISTING 



This section includes the BASIC computer code that was used to 
calculate the predicted heat fluxes at the exterior and interior surfaces of the 
roof with roof-spray cooling in Bryan, Texas. Although the code includes the 
data for many input parameters in the text, the code could be easily modified 
to read these parameters from an input data file or to prompt the user for the 
input values. The code could also be easily translated into FORTRAN if 
desired for compatibility with other software or data files. 

The following code calculates the predicted heat fluxes at the exterior 
and interior roof surfaces of a test building that Somasundaram et al. [3] used 
for roof-spray cooling experiments. In this example combined convective and 
radiative heat transfer coefficients were used, but the code allows the use of 
convective-only heat transfer coefficients if the radiative heat fluxes are not 
set equal to zero. The external convective heat transfer coefficient, h^, can 
also be varied if needed to account for large variations in wind speed, v. 

The iterative techniques used in this code proved to be effective in the 
required simulation. However, the speed of the solution could be improved 
through the application of more sophisticated bracketing and convergence 
techniques. 
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REM ROOF-SPRAY COOLED ROOF AT BRYAN, TEXAS (IN FILE 
ATMYBRTX BAS) 

REM HEAT FLOW CALCULATION BY TRANSFER FUNCTION 
METHOD 

OPEN ”C:ATMBRTX.PRN" FOR OUTPUT AS #1 

WRITE #1, ’’Hour", "OMEGA”, ”Tset”, ”Ta”, ”Twet”, ”TrO”, "Try”, 
"Qsolar", "Qevap", "QcondO", "Qe", "QoutEXP", "QinEXP", 
"qCONVo", "pWR", "pW”, ” " 

LET T=0 

5 LET H=T 

6 IF H>=24 THEN H=H-24 

IF H>=24 GOTO 6 

DIM Tr(-6 TO 120), TE(-6 TO 120), QE(-6 TO 120), Try(-6 TO 120), 
1(0 TO 120), QEXPout(0 TO 120), QEXPin(0 TO 120) 

REM Tdry-bulb, Twet-bulb, PRESSURE, HUMIDITY, AND SOLAR 
VALUES FROM TYPICAL METEOROLOGICAL YEAR 
WEATHER DATA AVERAGED FOR PERIOD 17-30 JULY FOR 
WACO, TEXAS 

REM WEATHER FILES OBTAINED FROM ACROSOFT 

INTERNATIONAL, INC., DENVER, CO AND PROCESSED BY 
MICRO-DOE2 SOFTWARE 

REM EXPERIMENTAL DATA FROM SOMASUNDARAM et al. [1988] pp. 
1103-1104 

IF H=0 THEN TW=69: TAMB=79.3: P=29.45: W=.031: QSOL=0; 
QEXPout=-3: QEXPin=.45 

IF H=1 THEN TW=69: TAMB=77.8; P=29.47; W=.0136; QSOL=0; 
QEXPout=-3: QEXPin=.34 

IF H=2 THEN TW=68.9: TAMB=76.5: P=29.47: W=.0138: QSOL=0: 
QEXPout=-3.15: QEXPin=.27 

IF H=3 THEN TW=68.7; TAMB=75.3: P=29.47; W=.014: QSOL=0: 
QEXPout=-3.I5: QEXPin=.24 

IF H=4 THEN TW=68.6: TAMB=74.6: P=29.48; W=.0I4: QSOL=0: 
QEXPout=-3.I5: QEXPin=.2I 

IF H=5 THEN TW=68.3; TAMB=73.7: P=29.48: W=.0I4; QSOL=0; 
QEXPout=-3.15: QEXPin=.17 

IF H=6 THEN TW=68.1: TAMB=72.9: P=29.49; W=.0I4: QSOL=0: 
QEXPout=-3.15: QEXPin=.14 

IF H=7 THEN TW=69.7: TAMB=76: P=29.51: W=.0145: QSOL=15.5: 
QEXPout=-2.47: QEXPin=-.21 

IF H=8 THEN TW=71.1: TAMB=79.1: P=29.51: W=.0I49: 
QSOL=104.9: QEXPout=-l . 1 : QEXPin=-.2 
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IF H=9 THEN TW=72.6; TAMB=82.2: P=29.51; W=.0153; 

QSOL=170.3: QEXPout=l .64: QEXPin=-.24 
IF H=10 THEN TW=72.7: TAMB=84 9: P=29.51: W=.0149; 

QSOL=202.6: QEXPout=3.83; QEXPin=-.24 
IF H=ll THEN TW=73.1: TAMB=87.4: P=29 51: W=.0146: 
QSOL=215.4: QEXPout=5.34: QEXPin=-.07 
IF H=12 THEN TW=73.3: TAMB=90.1; P=29.49: W=.0141: 
QSOL=218.2; QEXPout=6.59: QEXPin=.09 
IF H=13 THEN TW=73.1; TAMB=91.3: P=29.49: W=.0136: 
QSOL=213.1; QEXPout=7 12; QEXPin=.38 
IF H=14 THEN TW=72.9: TAMB=92.7; P=29.47: W=.0132; 

QSOL=206.9: QEXPout=5.62: QEXPin=.48 
IF H=15 THEN TW=72.9: TAMB=94; P=29.43; W=.0129; 

QSOL=205.4: QEXPout=5.2: QEXPin=.75 
IF H=16 THEN TW=72.5: TAMB=93.3: P=29.41; W=.0128; 

QSOL=199.1; QEXPout=3.7; QEXPin=1.09 
IF H=17 THEN TW=72.1: TAMB=92.7: P=29.41: W=.0126; 

QSOL=178.7: QEXPout=3.3: QEXPin=l 06 
IF H=18 THEN TW=71.5; TAMB=92: P=29.41: W=.0122: 

QSOL=124.2: QEXPout=l .64: QEXPin=1.06 
IF H=19 THEN TW=70.8; TAMB=88.7: P=29.41: W=.0124: 

QSOL=4.7: QEXPout=-. 13: QEXPin=1.35 
IF H=20 THEN TW=70.3: TAMB=85.9: P=29.42: W=.0127: QSOL=0: 
QEXPout=-l.l: QEXPin=1.16 

IF H=21 THEN TW=69.8: TAMB=82.8; P=29.43: W=.0129: QSOL=0: 
QEXPout=-2.33: QEXPin= 99 

IF H=22 THEN TW=69.5: TAMB=81.7; P=29.44: W=.0131: QSOL=0: 
QEXPout=-2.74: QEXPin=.75 

IF H=23 THEN TW=69.1: TAMB=80.4; P=29.45; W=.0130: QSOL=0: 
QEXPout=-2.86: QEXPin=.49 

IF H=24 THEN TW=69: TAMB=79.3: P=29.45: W=.031: QSOL=0: 
QEXPout=-3: QEXPin=.45 
REM END OF TMY DATA 

P=P*3386.4: REM CONVERT PRESSURE FROM INCHES HG TO Pa 
Pw=(P*W)/(.62198+W): REM CALCULATE PARTIAL PRESSURE OF 
WATER IN MOIST AMBIENT AIR 
LET QSOL=.6*QSOL*3.155: REM APPLY ROOF ABSORPTANCE 
AND CONVERT FROM Btuh/sf TO W/m**2 
Twet=(TW-32)/l. 8+273. 15: REM CONVERT WET BULB TEMP 
FROM F TO K 

TAMB=(TAMB-32)/1. 8+273. 15: REM CONVERT DRY BULB TEMP 
FROM F TO K 
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REM FROM ROOF GROUP #166 BUILD ON TRANE TRACE 600 

ROOF GENERATOR UTILITY BASED ON SOMASUNDARAM 
et al. TEST BUILDING 

REM DESCRIPTION OF ROOF, CONDUCTION TRANSFER 
FUNCTION COEFFICIENTS 

BO = 4.40399E-09*5.6783: REM CONVERT Bn AND Cn TO SI UNIT 
B1 = 1.461269E-05*5.6783; B2 = 2.489753E-04*5.6783 
B3 = 5.109677E-04*5.6783; B4 = 2.08079E-04*5.6783 
B5 = 1.789627E-05*5.6783: B6 = 2.828292E-07*5.6783 
Cn = .001*5.6783: DO = 1 

D1 = -1.974663; D2 = 1.273348: D3 = -.3056189 
D4 = .02383278: D5 = -4.3 12459E-04. D6 = 1.107164E-06 
REM TIME LAG = 9 HOURS, Utable = .061*5.6783, W/(m**2*K) 
REM DF = .227, DECREMENT FACTOR, AMPLITUDE REDUCTION 
FACTOR FROM OUTSIDE HEAT FLUX TO INTERIOR LOAD 
Ho = 22.7: REM EXT CONVECTIVE AND RADIATIVE HEAT 
TRANSFER COEFFICIENT, W/(m**2*K) 

Hi = 8.35: REM INT CONVECTIVE AND RADIATIVE HEAT 
TRANSFER COEFFICIENT, W/(m**2*K) 

EPS = .85; REM EMISSIVITY OF ROOF 
SIG = 5.67E-08; REM STEFAN-BOLTZMAN CONSTANT, 
W/(m**2*K**4) 

Rr = 2.079: REM THERMAL RESISTANCE OF ROOF, m**2*K/W 
TSET = 299.81: REM AMBIENT TEMPERATURE OF INSIDE AIR 
PLENUM, K 

LET Try(-l) = 280; REM PROVIDES STARTING TRIAL TEMP FOR 
HOUR 0 

IF T>=24 THEN Try(T)=Try(T-24)- 1.0001; GOTO 20 
Try(T) = Try(T-l) - 20 

20 Fry = Try(T)*l 8: CONVERT TEMP FROM K TO F 

REM CALCULATE SATURATION PRESSURE OF WATER ON 
ROOF SURFACE FROM MATHUR [1988] 

PARTO - 670.2012: PARTI = -5.325521 *Fry 

PART2 = 0159464*Fry**2: PART3 = -.00002134061 *Fry**3 

PART4 = 1.077853E-08*Fry**4 

Pwr = (PART0+PART1+PART2+PART3+PART4)*6894.8: REM Pa 
QEVAP = 073814*(Pwr-Pw): REM CALCULATE EVAPORATIVE 
HEAT TRANSFER DUE TO PRESSURE DIFFERENCE OF 
WATER VAPOR AT ROOF SURFACE AND IN AMBIENT AIR 
REM Tr(T) = TAMB+Rr*(-QEVAP-QRAD-QCONV+QSOL) 

QRADo = EPS*SIG*(Try(T)**4-(TAMB-3)**4) 

QCONVo = Ho*(Try(T)-TAMB) 
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Tri = Try(T)+Rr*(QRADo+QEVAP+QCONVo-QSOL); REM CALC 
INT ROOF TEMP, K 

QCONDo = -QRADO-QEVAP.QCONVo+QSOL 
QCONVi = Hi*(Tri-TSET) 

QRADi = EPS*SIG*(Tri*M-TSET**4) 

REM QCONDi = QCONVi=QRADi; REM INT SURFACE ENERGY 
BALANCE 

Tr(T) = Tri+Rr*(QCONVi+QRADi) 

IF ABS(Tr(T)-Try(T))< l GOTO 40: REM CHECK FOR 
CONVERGENCE OF TRIAL ROOF TEMP WITH 
CALCULATED TEMPERATURE 

REM SET STEP SIZE TO CALCULATE OUTSIDE ROOF TEMP 
IF Tr(T)>285 AND Tri>0 THEN I(T)=.0001: GOTO 30 
IF Tr(T)>270 AND Tri>0 THEN I(T)=.001: GOTO 30 
IF Tr(T)>0 AND Tri>0 THEN I(T)=.01; GOTO 30 
IF Tr(T)<0 OR Tri<0 THEN I(T)=.l 

30 LET Try(T) = Try(T)+I(T): REM ADD ITERATIVE STEP TO 
PREVIOUS TRIAL ROOF TEMPERATURE 
IF Try(T)>Tr(T-l)+20 THEN PRINT "ERROR ON ITERATIVE 
LOOP": GOTO 50 
GOTO 20 

40 REM SOL-AIR TEMPERATURE WILL BE TAKEN AS 

CALCULATED ROOF SURFACE TEMPERATURE 
LET TE(T) = Tr(T)-273.15: REM CONVERT ROOF TEMP FROM K 
TO C 

REM CALCULATE PREDICTED LOAD ON INTERIOR FROM ROOF 
BY TRANSFER FUNCTION METHOD 
QEl=B0*TE(T)+Bl*TE(T-l)+B2*TE(T-2)+B3*TE(T-3) 
QE3=B4*TE(T-4)+B5*TE(T-5)+B6*TE(T-6) 
QE3=Dl*QE(T-l)+D2*QE(T-2)+D3*QE(T-3) 
QE4=D4*QE(T-4)+D5*QE(T-5)+D6*QE(T-6)+(TSET-273.15)*Cn 
QE(T) = QE1+QE2-QE3-QE4 
REM WRITE OUTPUT DATA TO FILE 

WRITE #1, T, W, TSET, TAMB, TWET, Tr(T), QSOL, QEVAP, 

QCONDo, QE(T), QEXPout(T), QEXPin(T), QCONVo, Pwr, Pw 
LET T = T+1: REM START NEXT HOUR SIMULATION 
IF T>120.5 GOTO 50: REM STOP SIMULATION AFTER 120 HOURS 
GOTO 5 

50 REM END OF WHILE LOOP 

CLOSE #1: REM CLOSE OUTPUT FILE ATMYBRTX.PRN 
END 
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REM ******* * * NOMENCLATURE* * *************************** 
Bn = CONDUCTION TRANSFER FUNCTION COEFFICIENTS 

Cn = CONDUCTION TRANSFER FUNCTION COEFFICIENTS 

Dn = CONDUCTION TRANSFER FUNCTION COEFFICIENTS 

EPS =EMISSIVITY 

Fry = TRIAL ROOF TEMPERATURE, DEGREES R 

H = HOUR OF DAY, MODULUS 24 CLOCK, HOURS 

Hi = INTERNAL CONVECTIVE(/ RADIATIVE) HEAT TRANSFER 

COEFFICIENT, W/(M**2*K) 

Ho = EXTERNAL CONVECTIVE(/ RADIATIVE) HEAT TRANSFER 
COEFFICIENT, W/(M**2*K) 

P = TOTAL PRESSURE OF MOIST AIR, INCHES HG, PSIA OR Pa 

Pw = PARTIAL PRESSURE OF WATER IN MOIST AIR, Pa 

Pwr = PARTIAL PRESSURE OF WATER AT ROOF SURFACE (Psat 

AT Troof), Pa 

Rr = THERMAL RESISTANCE OF ROOF, m**2*K/W 

QCONDo - CONDUCTIVE HEAT FLUX INTO ROOF FROM OUTER 
SURFACE, W/m**2 

QCONVi = CONVECTIVE HEAT FLUX AT BOTTOM OF ROOF TO 
ROOM, W/m**2 

QCONVo = CONVECTIVE HEAT FLUX AT TOP ROOF TO AMBIENT, 
W/m**2 

QE(T) = PREDICTED HEAT FLUX INTO ROOM, W/m**2 
QEVAP = EVAPORATIVE HEAT FLUX FROM MOIST ROOF, W/m**2 
QRADi = RADIATIVE HEAT TRANSFER FROM ROOF INNER 
SURFACE TO AMBIENT ROOM, W/m**2 
QRADo = RADIATIVE HEAT TRANSFER FROM ROOF OUTER 
SURFACE TO AMBIENT, W/m**2 

QSOL = SOLAR RADIATION, INCIDENT AND ABSORBED, W/m**2 
SIG = SIGMA-STEFAN-BOLTZMAN CONSTANT, 5.6697E-08 
W/(m**2*K**4) 

T = NUMBER OF HOURS STARTING AT MIDNIGHT ON FIRST 

SIMULATION DAY, H 

TAMB = AMBIENT DRY BULB TEMPERATURE, F, C OR K 
TE(T) = CALCULATED EQUIVALENT ROOF TEMPERATURE =(Tr), C 
Tr(T) = ROOF SURFACE TEMPERATURE AT HOUR T, C OR K 
Tri = CALCULATED INSIDE ROOF SURFACE TEMPERATURE, K 
Try(T) = TRIAL ROOF SURFACE TEMPERATURE, K 
TSET = AMBIENT ROOM/ AIR PLENUM TEMPERATURE, C OR K 
TW = AMBIENT WET BULB TEMPERATURE, F, C OR K 
W = OMEGA-HUMIDITY RATIO, MASS OF WATER PER UNIT 

MASS OF DRY AIR 



APPENDIX B 

UNCERTAINTY ANALYSIS 



To estimate the accuracy of the predicted results it is necessary to 
quantify the uncertainty of the individual variables and its effect on the 
uncertainty of the model results. In this thesis the result is the prediction of 
the heat gain through the roof, q^. The variables include ambient temperature, 
humidity ratio, solar heat flux at the roof surface, inside and outside 
convective heat transfer coefficients, thermal resistance of the roof, and net 
radiative heat exchange with the sky. 

The uncertainty analysis method chosen is the sequential perturbation of 
the data reduction program described in Appendix A. This method is 
described in detail by Moffat [19]. This method is based on determining the 
effect of each variable's uncertainty on the baseline predicted result. After the 
result of the baseline data set are calculated results are determined "once more 
for each variable, with the value of the variable increased by its uncertainty 
interval (and all other [variables] returned to their baseline values)". [19] The 
difference between the returned perturbed results and the baseline represents 
the contribution of each variable's uncertainty on the uncertainty of the result's 
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overall uncertainty The uncertainty of the result is determined by squaring 
each contribution, summing, and taking the square root. 

The expression used to calculate the uncertainty W of a function is 



where Xj is any of the variables which quantify the function. 

The uncertainty interval for each variable has been estimated based on 
Somasundaram et al. [5] and Houghten et al.[l] . The results of the sequential 
perturbation for the roof-spray cooled concrete roof at Pittsburgh are 
provided in Table B.l. 

Table B.l Uncertainty Analysis by Sequential Perturbation for Concrete 




(B-1) 



Roof at Pittsburgh 



Case/ 

Perturbed 

Variable 



Uncertainty 

Interval 



Perturbed Difference From 
Result Baseline 



(W/m') (W/m') 



Baseline 



-4.73 

-4.28 

-4.72 

-4.27 

-4.71 

-4.28 

-4.42 

-5.42 

-5.13 



0.5 K 



0.45 

0.01 

0.46 

0.02 

0.45 

0.31 

0.69 

0.39 



0.001 



q„u. and “ 

hi 

ho 

R, 



10 W/m^ 

2 W/(m" K) 
4 W/(m^ K) 



0.05 KAV 




0.5 K 
10 W/m^ 
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The calculation of the total uncertainty, W, is shown by 
W= {(0.45)'+(0.01)'+(0.46)'+(0.02)'+(0.45)''+(0.31)'+(0.69)'+(0.39)'}®' 
which yields an uncertainty, W, of 1.2 W/m^. Similar calculations for the 
other roofs yielded uncertainties of 0.2 W/m^ for the predicted 0.40 W/m^ 
predicted average roof heat gain for the spray cooled roof at Bryan, Texas, 1.2 
W/m^ for the predicted 2.4 W/m^ predicted average for the dry roof at Bryan, 
and 1.4 W/m^ for the predicted -4.2 W/m^ average for the spray cooled pine 
plank roof at Pittsburgh. 
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